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On the basis of the Hubbard model, we present the formulation of antiferromagnetism in electron- 
doped cuprates using the fluctuation-exchange approach. Taking into account the spin fluctuations 
in combination with the impurity scattering effect due to the randomly distributed dopant-atoms, we 
investigate the magnetic properties of the system. It is shown that the antiferromagnetic transition 
temperature, the onset temperature of the pseudogap formation, the single particle spectral density, 
and the staggered magnetization obtained by the present approach are in very good agreement with 
the experimental results. The distribution function in momentum space at very low temperature is 
observed to differ significantly from that of the Fermi liquid. Also, we find zero-energy peak in the 
density of states (DOS) of the antiferromagnetic phase. This DOS peak is sharp in the low doping 
regime, and disappears near the optimal doping where the AF order becomes weak. 
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I. INTRODUCTION 

The electron-doped cuprate high-temperature super- 
conductors (EDCHTS) such as Nd^Ce^CuCU (NCCO) 
and Pr 2 _ 2; Ce 2 ;Cu04 (PCCO) are in the antiferro- 
magnetic (AF) state within a wide doping range at 
low temperature ^2^4 The angle-resolved photoemission 
spectroscopy (ARPES) measurements^!^ have revealed 
the Fermi surface (FS) evolution with doping, provid- 
ing some microscopic information about the antiferro- 
magnetism in these compounds. On theoretical as- 
pect, much work has been carried out to understand 
the electronic properties of the EDCHTS by using cither 
Hubbard 8 i 9 i 10 i n i 12 i 13 i 14 or t - J type models^ But there 
are difficulties with the existing theories in describing 
the AF phase transition and the FS evolution with dop- 
ing in a consistent manner. To explain the FS evolution 
and the staggered magnetization as a function of doping, 
some authors have to assume a doping-depcndcnt Hub- 
bard J7i2iii Though the mean-held theory (MFT) yields 
qualitative explanation for the FS evolution, the calcu- 
lated AF transition temperature Jjv is well known to be 
too highiiiii^ In order to reduce the magnitude of Tjv, the 
self-consistent random-phase approximation which takes 
into account the effect of spin fluctuations has been em- 
ployed, but the obtained Tm is a non-monotonic function 
with increasing doping concentration rrjii contrary to the 
experimental observation. Therefore, a quantitative un- 
derstanding of the AF phase transition in EDCHTS re- 
mains a challenging issue. 

The high-temperature superconductors (HTS) based 
on cuprates are considered as typical quasi-two- 
dimensional (Q2D) electron systems. At the early stage 
after the discovery of the HTS, the antiferromagnetism 
in such systems has been studied by the random-phase 
approximation (RPA) with the Hubbard model piSiiLi&iS 
In this approach the spin fluctuations are treated as per- 
turbation over the mean-held AF state. For the Q2D sys- 
tems, however, this treatment for the AF state at finite 



temperature may not be adequate because the number 
density of the long-wavelength fluctuating spins diverges 
logarithmically when the dimensionality approaches the 
limit 2. It is therefore desirable to treat the AF order 
and the spin fluctuations on an equal- footing manner. 

Very recently, two of us have extended the fluctuation- 
exchange (FLEX) approach based on the Hubbard 
model, which was designed for studying the spin fluc- 
tuations in a two-dimensional strongly-correlated elec- 
tron system^ to investigate the antiferromagnetism be- 
low Tat in NCCO or PCCO^i In the formalism, the 
single-particle Green's function, the AF order parame- 
ter, and the charge and spin response functions are all 
treated self-consistently. We found that, in order to 
quantitatively account for the AF properties, not only 
the spin fluctuations, but also the impurity scattering ef- 
fect due to the randomly distributed dopant-atoms are 
necessary to be taken into consideration. As we know, 
while the Ce atoms are doped to modulate the carrier 
density they should simultaneously introduce disorder 
which influences the property of electrons in CuO lay- 
ers. By incorporating the impurity scattering into the 
FLEX formalism, we have shown, with the Hubbard in- 
teraction U = 8, that the numerically obtained Tjy, the 
onset temperature of pseudogap due to spin fluctuations, 
the single particle spectral density near FS, and the stag- 
gered magnetization m in the AF phase as a function of x 
are all in good agreement with the relevant experimental 
measurements. 

In this paper, we want to present the detailed FLEX 
formulation of the AF theory in the EDCHTS which is 
omitted in Ref. l2ll and show some improved and new 
results as well. By improving the treatment of interlayer 
coupling and slightly adjusting the parameters, better 
agreements with experiments can be achieved. As a new 
result, the single-particle distribution function has been 
obtained. Its behavior is found to be significantly dif- 
ferent from that of the Fermi liquid at low temperature. 
In addition, the zero-energy peak is found in the density 
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of states (DOS) of the AF phase. This peak is sharp in 
very low doping regime. The peak becomes suppressed 
as the doping concentration gets larger and completely 
disappears near the optimal doping. 



II. FORMALISM 

We start with the Hubbard model including the impu- 
rity field: 
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where £g = — \l with as the tight-binding energy 
spectrum, /i the chemical potential, ct (cg CT ) is the elec- 
tron creation (annihilation) operator with spm-cr(= ±f 
for up and down spins, respectively) and momentum k, 
n a a = y~!r ct Ct* i — , and N is the total number of the lat- 
tice sites. The quantity Vi (assumed to be a constant) in 
the last term represents the potential acting on the elec- 
tron with density ni at the impurity site /. In NCCO 
or PCCO, the concentration Ni of the dopant Ce atoms 
is the same as the electron-doping concentration x. The 
electron system under consideration is a Q2D one with 
weak interlayer coupling. 

Our main task here is to give a scheme for determining 
the electronic Green's function G. In the normal state, 
G is related to the self-energy E via 



G(k,iuj n ) 
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where uj n is the fermionic Matsubara frequency. For 
brevity, we hereafter will use k — (k,iu> n ) for the ar- 
gument unless other illustration for special case. The 
formalism for G is to find out the self-energy X as a 
functional of G. For constructing the self-energy dia- 
grams, it is the usual way to separate the [/-interaction 
term in Eq. Q into channels of charge density and spin 
fluctuations^ 
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where riff = ?/„■ 



<T\ 



+ riff 



52k c l+tf c kl ( S q = S P"> are the char S e density, 
longitudinal-spin and transverse-spin fluctuation opera- 
tors, respectively. Then, diagramming of the self-energy 
under the FLEX approximation is to piece together these 
interaction vertices in the right-hand-side (r.h.s.) of Eq. 
(PI in each channel. 

Having obtained Eq. J3J), we also need to take into ac- 
count the effect due to the interlayer spin coupling. For 
NCCO or PCCO in which the CuO layers are stacked in 
a staggered way as shown in Fig. 1, the AF spin coupling 
between the nearest CuO layers is highly frustrated, and 
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FIG. 1: Crystal structure of NcbCuGi. The arrow shows the 
spin orientation on the Cu site. 



is zero in case of ideal stacking. The local lattice dis- 
tortions or any other asymmetry may result in a weak 
AF spin coupling between the nearest layers. In the case 
of ideal stacking, the coupling between the next-nearest 
CuO layers is then dominant. We here suppose the ideal 
stacking and introduce the next-nearest-layer spin cou- 
pling Hamiltonian, 



(4) 



where i^'-summation is over the next-nearest CuO layers. 
By adding this z-direction coupling to the Hubbard-C/ 
term, the interaction is then given by 

H ^ = ^E^-^-^E^^^v^^]' ( 5 ) 



with 



V(q) = U - J z cosq z . 



(6) 



The argument q z implies that the length unit in z- 
direction is the next-nearest interlayer distance. The 
weak interlayer coupling means \J Z \ -C U. Having the 
constant J z to mimic the interlayer AF coupling, we then 
neglect the interlayer electron hopping. The energy 
in Eq. 0]) is hereafter considered as a function of two- 
dimensional momentum k. 

In the following subsections A and B, we construct 
the self-energy diagrams in the normal and AF states, 
respectively. 



Normal state 



In the normal state, following the FLEX 
approximation, 2 4 we can easily get the self-energy 
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FIG. 3: Self-energy due to the spin fluctuations in the normal 
state. The point-wavy line represents the interaction —V(q). 
(a) Ring diagrams representing the longitudinal-spin fluctua- 
tions, and (b) transverse-spin fluctuations. 



FIG. 2: Diagrams for self-energy due to the charge-density 
fluctuations in normal state. The wavy line represents the on- 
site Hubbard U. The dashed single line with a cross denotes 
the impurity scattering, (a) Lowest order (Hartree) terms, 
where a double counting (d.c.) without the impurity scat- 
tering should be subtracted, (b) ring diagrams representing 
the charge-density fluctuations, (c) vertex corrections due to 
the impurity scattering for the particle-hole and (d) particle- 
particle channels. 



diagrams in the presence of impurities. Firstly, we 
consider the contribution from charge-density fluctua- 
tions represented by the first term in r.h.s. of Eq. JSJ. 
Its contribution to the self-energy is diagrammatically 
shown in Fig. [3 The diagrams in Fig. |2Ia) come from 
the Hartree term. The first one includes the diffuson 
correction due to the particle-hole propagator under 
the impurity scattering. The second one contains the 
Cooperon correction. The contribution to the self-energy 
in the absence of impurity scattering is subtracted be- 
cause of a double counting and the remaining constant 
absorbed in the chemical potential. The ring diagrams 
in Fig. Elb) represent the charge-density fluctuations 
with the decoration of the impurity scattering. The 
bubble in Fig. |2Ib) is the charge susceptibility defined 
by 

x(q,t - r') = --^(T T n^{T)n-^(T')) bub (7) 

where the subscript bub means the bubble diagram with 
impurity insertions. The impurity insertions or vertex 
corrections in the above diagrams are illustrated in Figs. 
Wi c ) and 1(d). They can be expressed by the same for- 
mula, 



T(q,iuj n ) = [1- N t vfX(q,iuj n )}- 



(8) 



with 



X(q,iu n ) = G(k, iu) n )G(k + q, iu n + iQ, m ), (9) 



where q = (q,ifl m ) with f2 m the bosonic Matsubara fre- 
quency. With the vertex correction, in the frequency- 



space, the charge susceptibility \ is given by 
X(v) = jj^2 T (q,iu n )X(q,iuj n ). 



(10) 



The part of the self-energy due to charge-density fluctu- 
ation is then expressed as 
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is the polarization potential due to the charge-density 
fluctuations. 

Analogously, we can treat the spin fluctuations given 
by the last term in r.h.s. of Eq. (JSJ). Figures 3(a) and 
3(b) are the self-energy diagrams due to longitudinal- 
spin and transverse-spin fluctuations, respectively. These 
diagrams are almost the same as Fig. |2tb) except the 
bubbles with different meanings. In Fig. 01 a), the bubble 
is the longitudinal-spin susceptibility defined by 



X z (q,r - t') ee --(T T S%(T)S%~(T')) bubl 



(13) 



while in Fig. Ofb), it is the transverse-spin susceptibility 



xH^t-t') ee --(T T S ? (r)S+(r')) bu 



(14) 



It is easy to see that all three susceptibilities defined 
above are the same in the normal state. Therefore, we 
hereafter simply use the notation x(<?) for them without a 
superscription. The contribution to the self-energy from 
the spin fluctuations can then be expressed as 

S s (fc) = ~ G ( k + l)il VS (l) U 2 X (q)}T 2 (q, ico n ), 
i 

(15) 
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FIG. 4: Self-energy due to the impurity scattering. 



where 



V S (q) 



1 + V(q) X (q) 



(16) 



is the potential produced by the spin polarization. The 
term — U 2 x{q) in Eq. H15fl represents the subtraction of 
a double counting in second order of U. 

The process of the single electron scattering off impuri- 
ties is shown in Fig. Its contribution to the self-energy 
is given by 




FIG. 5: (Color online) Static spin susceptibility x(<7, 0) x (U + 
| J, |) at T = T N and x = 0.1. 
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The total self-energy is then given by the summation, 
£ = £ c + E s + £\ Note that by setting v l = and 
J z = 0, £ reduces to the previous result of the FLEX 
approximation^^ 

At temperature above but close to T/v, the most impor- 
tant contribution comes from the spin fluctuations repre- 
sented by Figs. 02a) and|3fb) since they give rise to the 
AF instability at T/v. The AF instability is character- 
ized by the singularity in V s (q). At q — (Q,0) = Q, the 
magnitude of x(o) reaches its maximum. At T = T/v, we 
have 1 + V(Q)x{Q) = 0. 

We here emphasize the importance of the interlayer 
coupling. Without interlayer spin coupling, i.e., J z = 0, 
the integral in the r.h.s. of Eq. (|15fl will diverge logarith- 
mically at q = Q and finite temperature, which means 
that no AF transition occurs at finite temperature for 
two-dimensional (2D) systems. Though J z is small, it 
controls the overall magnitude of the AF transition tem- 
perature T/v. 

The interlayer spin coupling becomes significant only 
at the singular point of V s (q). The Green's function 
G and the susceptibility x verv weakly depend on the 
weak interlayer coupling. By neglecting the weak q z - 
dependence in G and %, we can then integrate out 
Eq. (|15|) over q z immediately, and reduce the three- 
dimensional problem to a 2D one. Since the delicate 
q z -dependence of V s {q) comes from the denominator in 
Eq. I |lti|l . we neglect the small J z in V 2 {q) in the nu- 
merator. By integrating over q z , V s (q) is replaced by 



V s (q) 



U 2 x(q) 
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This result is just the average of V s (q) over — it < q z < ir. 
In our numerical calculation, \J Z \/U < 2%, the above 
approximation should be reasonable. 



Figure[S]shows a numerical result for x(Qi 0) at T = T/v 
for doping concentration at x — 0.1, and it indicates that 
the maximum of |x(<f, 0)| appears at {Qx,Qy) = (tTjTt)- 
Our numerical calculation also shows that such behavior 
for x(q> 0) remains unchanged in the whole doping range 
investigated here. This is different from the hole-doped 
case where the maximum appears at some incommen- 
surate wavevectors close to (7r,7r). The reason is that 
the Fermi surface can be nested at (Q x ,Q y ) = (tt, tt) for 
the electron doped sample within a large doping range. 
For the hole doped sample, the Fermi surface cannot be 
nested at {Q x , Qy) = (tt, tt) unless in the very low doping 
regime. As for the z-component, Q z = or ir depends on 
the sign of J z . Under our approximation where V s (q) is 
replaced by V s (q), the final result is independent of Q z . 
In NCCO and PCCO, Q z = is observed, which means 
a negative J z or ferromagnetic interlayer coupling. 



B. AF phase 

Below T/v, the system is in the AF state that is char- 
acterized by the staggered magnetization m. Now, the 
operator — > —Nm is a macroscopic quantity. The 

term of q = Q in the spin interaction in Eq. (J5J should 
be singled out, and be treated as (up to a constant) 
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where A = V(Q)m/2 is the order parameter. Though 
Eq. i|19[l looks like the mean-field result, we here also 
take into account simultaneously the fluctuation terms 
q Q in the interaction. Because of Eq. 1191) . we need 
to redefine the Green's function as a 2 x 2 matrix 



G(k, 



') = -<^>)4>')> (20) 
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FIG. 6: Off-diagonal self-energy E12. 



with ibt = (cl , act - ). Since the potential (±A) in 

Eq. (|19|l for opposite spin electrons has opposite sign, G 
so defined is independent of the spin. In the Matsubara 
frequency space, in terms of the self-energy E(fc) which 
is accordingly a 2 x 2 matrix, G(k) is written as 



G(k) = [iw„ - £ (k) - &{k)a 3 - E(fc)]" 



(21) 



where £ (3)( fc ) = (% ± ct 3 is the Pauli matrix, 

and the scalar quantity ito„ — £o(&) is actually multiplied 
by the unit matrix <jq that is dropped out for brevity. 
The four elements of G have the following relationship: 

G 22 (fc) = Gu(k + Q) 

G 2 i(fc) = Gia(fc) = G 12 (fe + Q). (22) 

Similarly, the elements of E have the following relation- 
ship: 



E 22 (fc) - E n (fc + Q) 

E a i(fc) - Ei 2 (fc) = Ei 2 (fc + Q). 



(23) 



By definition, the order parameter A can be written as 



A = 



v(Q) 

Nfi 



E G ^( fc ) 



(24) 



where (3 is the inverse of the temperature T. From Eq. 
()19J) with the impurity influence, the off-diagonal self- 
energy Ei 2 (fc) is obtained as in Fig. 6. 

For the diagonal self-energy, we need to extend the re- 
sults shown in Figs. 12141 for the normal state to that of 
the AF state. Since the diagonal and off-diagonal parts 
of the self-energy have a kind of harmonious relation- 
ship, they cannot be treated independently. Here, two 
points need to be noted: (1) Except the diagram shown 
in Fig. El are there any others that should be included 
in Ei 2 (fc) under the ring approximation for the diagonal 
part? The transverse-spin fluctuations corresponding to 
Fig. E{b) should include the collective modes, especially 
the Goldstone mode (GM) that is a consequence of the 
broken symmetry. The equation determining the GM 
should be consistent with the gap equation (J52J. This 
consistency gives rise to a constraint on the form of E 12 . 
We will show that Fig. EJis the only graph for E12 within 
the ring approximation for the diagonal self-energy. (2) 



What should be the form of the impurity insertion? The 
insertion of impurity scatterings in the transverse-spin 
fluctuation is between the opposite spin electrons, while 
in Fig. [5] it is between the same spin electrons. In or- 
der to ensure the consistency between GM and Eq. 124fl , 
these insertions should be the same, and the scattering 
processes including Gi 2 will be excluded. In the follow- 
ing, we give detailed explanations to these two points 
with and without impurity scatterings. 

1. Without impurity scattering (i.e., Vi = 0). In the 
AF state, since the transverse-spin fluctuations of mo- 
menta q and q + Q are correlated, we need to consider a 
matrix x for the transverse-spin susceptibilities. For the 
purpose, we denote the transverse-spin operators by 



(25) 



where = (for fi = 1), or (J(for /i = 2). The matrix 
elements are then defined by 

XU(Q,T- r') = -±(T T S^(q,T)St(q;T')) bub . (26) 
After Fourier transform from r to i£l m , they are given by 
Xilfa) = ^E[ G "( fc ) G "( fc + 9) -G 12 (k)G 12 (k + q)\ 

k 

X*a(9) = m E[ G "( fc ) G i2(fc + q) ~ G 12 (fc)G n (fc + q)] 

k 

and X22O) = Xu(q + Q), xliio) = Xizio)- In terms 
of X t (q)i we can then express the transverse-spin polar- 
ization potential V t (q) [stemming from a similar bubble 
summation as in Fig. 12b] as 

V\q) - [1 + ViqWiqT'nqWiqWiq), (27) 

which is a 2 x 2 matrix too. Note that V(q) is the scalar 
V(q) as given by Eq. @) if Q z = 0, or a matrix V(q) — 
U — cr^Jz cos q z if Q z = 7r. Here, an important point is 
the divergence of V t {q) at q = Q [where x\i{Q) = 0]j 
that is 



l + V n (Q) X t 11 (Q) = 



(28) 



which means the existence of the Goldstone mode. To see 
the constraint on the off-diagonal diagrams, we express 
Eq. (|28|l in the following form. Since Vu(q) = V(q), we 
drop the subscripts in the following. Using the relation- 
ship between G and E, we have 



V(Q)x\i(Q) 



V(Q) 

N/3 

v(Q) 

N/3 



EGi 2 (fc)/E 12 (fe) 



(29) 
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FIG. 7: Contribution to the off-diagonal self-energy £12 in 
RPA. Under the present approximation, this term should be 
excluded. 




k t k+Q 



FIG. 8: A form for the impurity insertions in the off-diagonal 
self-energy £12. Under the present approximation, this vertex 
correction is not consistent with the Goldstone mode, and 
should be excluded. 



By comparing the last equation above with the gap equa- 
tion one gets 

Ei 2 (fc) = A. (30) 

In other words, if one deals with £12 (A;) improperly, Eq. 
(|28|l cannot be satisfied. Therefore, under the ring ap- 
proximation for the diagonal self-energy, A is the only 
term for £12 (k). 

It is suggestive to compare the above result with that 
from RPA. Within RPA the contribution from the cou- 
pling between the transverse spin fluctuations Si(q) and 
S\(q) is taken into account in £12 (k). Such a diagram 
is shown in Fig. However, this contribution should 
be excluded under the present approximation. Accord- 
ingly, the similar contribution from such coupling to the 
diagonal part should be excluded as well. The diago- 
nal self-energy due to the transverse-spin fluctuations is 
given by 

st ( fc ) = -jf'EW + tiWxib) - tfVute)], ( 31 ) 

where again the term — U 2 x\i(q) represents the subtrac- 
tion of the double counting in the second order of U. 

2. In the presence of impurity scatterings. The off- 
diagonal self-energy £12 represented by the diagram in 
Fig. [21 is written as, 

£ 12 (fc) = Ar(Q,iw n ). (32) 

with r(Q, iu! n ) as vertex correction due to impurity scat- 
terings. Note that the spins on the two sidelines of the 
vertex T(Q, iui n ) in Fig. EJare the same. In the AF state, 
besides the Green's function Gu, one may consider the 
insertions including G12 in T(Q,iuj n ). If G12 exists in 
the vertex correction in Fig. El it must appear in both 
sides of an impurity line as a couple. Since their spins 
are the same, one obtains a term of G\2{k)Gi2{k) in the 
fc-summation in X(Q,iuj n ). In the following, however, 
we show that these terms including Gi 2 should be ex- 
cluded in order to maintain the consistence between the 
Eq. (|32|l and the Goldstone mode. Analogous to what is 
done above, we consider the transverse-spin polarization 
potential V t (q) in the presence of the impurities. The 
potential ^(q) is given by Eq. I|27|l but with the vertex 
correction to the susceptibility included. We here only 



consider the case of q = Q, where V t (q) is singular. It 
means 

V(Q) X \i(Q) = ^Y / T 11 (Q,tLu n )[G 11 (k)G 22 (k)-G 2 12 (k)] 

= ^E G i2(fc)r u (Q,io; n )/£ 12 (A ; ) 
= -1, 

where Fn(Q, iuj n ) is the impurity insertion in the suscep- 
tibility Xii(Q)- Comparing this result with Eq. (|24p. we 
have £i2(fc) = AFn(Q, icu n ). On the other hand, from 
Eq. il3*2"|) . we get Tu(Q,iuj n ) — Y{Q,iuj n ). However, 
Tn{Q,iL0 n ) is the impurity insertion with the opposite 
spins in the two legs of the ladder. If there is G12, one 
should get a term — G\2{k)G\2{k) in T\\{Q, iuj n ), whose 
sign is different from that in T(Q,iuj n ). Therefore, those 
terms including G12 should be excluded from the vertex 
corrections. 

The approximation for the off-diagonal self-energy is 
very constrained by the consistence between the gap 
equation and the Goldstone mode. For example, shown 
in Fig. [HI is a possible case for the impurity insertions in 
£12. But this kind of vertex correction is not consistent 
with the Goldstone mode under the present approxima- 
tion, and therefore should be excluded. 

With these two points clarified, we can then give the 
expression for the self-energy in the AF state. We start 
with the expression for the vertex corrections. As dis- 
cussed above, the functional form of the vertex correction 
is given by Eq. I|27[) but with the Green's function G un- 
derstood to be Gu, rfo Wn) = [1 — N i v 2 X(q,iuj n )]^ 1 , 
with 

Xiq^^^^^Gn^Gnik + q). 

k 

With the vertex correction, the longitudinal \ z an d 
transverse x spin susceptibilities are given by 

X Z (q) = -rj^2 r (<l' iuj n)X + (q,iuj n ), (33) 

n 

X*i(?) = ^r(q,iw n )X-(q,iw n ), (34) 

n 

Xi2(q) = h^^{^ iuj n)X {q,iu n )T{q + Q,iw n ), (35) 
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with 



X±(q,iu n ) = - ^[Gn(fc)Gn(fc + q) ± G 12 (k)G 12 (k + q)\ 



X (q, iij n ) = i ^[G n (fc)G 12 (fc + 9 ) - G 12 (fc)G n (fc + g)]. 



The expression for S u reads 

Su(*) - -— ^G 11 (fc + 9 )[r 2 ( <z , l ^„)-l] + E l ^„) 

-^«E Gll ( fc + 9)^ff(<z)r 2 (g, iw n ), (36) 

with K eff (g) = V c (g)/2 + V'(q)/2 + V&(g) - tfV u fa) 
and 



£/V(g) 
i-£/x 2 (<zr 
^ 2 (g)x z (g) 



(37) 
(38) 
(39) 



The term V c (q) is due to the charge-density fluctuations, 
while the y z (g)-term stems from the longitudinal-spin 
fluctuations. At low T with finite A, since V^ 1 (q) di- 
verges at q = Q, the most important contribution to 
the diagonal self-energy comes from the transverse-spin 
fluctuations. The corresponding diagram is an extension 
of Fig. I3b). The other elements of £(fc) are given by 
£ 22 (fc) = En(fc + Q), and £ 21 (fc) = Eia(Jfe). 

Here again, we discuss the interlayer coupling J z that 
appears in V z (q) and V*i(q). For V z (q), since its sen- 
sitive dependence on J z comes from the denominator 
that vanishes at T/v in r.h.s. of Eq. (|38l) . we ne- 
glect the J z -dependence in the numerator and average 
it over — n < q z < tt as what have been done for nor- 
mal state. By the same consideration, we can deal with 
Vf\ (q) . But the situation for (q) is a little more com- 
plicated, we consider it in Appendix. In our previous 
calculation,— the interlayer spin coupling was taken into 
account artificially by introducing a weak q 2 -dependence 
in these denominators, which corresponds to the expan- 
sion of cosq z w 1 — g 2 /2 in the denominators to order 
of q z but neglecting the J z -dependent constant term. 
The physics is that the long- wavelength spin-fluctuations 
should be treated as three-dimensional, which is the same 
as present consideration. 



III. NUMERICAL RESULTS 

In the present calculation, we take the nearest- 
neighbor (NN) hopping t as the energy unit (i.e., t = 1) 
and chooser' = -0.25 (next NN hopping), t" = 0.1 (third 
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FIG. 9: (Color online) Tn as a function of doping concentra- 
tion x. Solid (vi/t — 1.55) and dashed (w, = 0) curves are the 
present calculated results. The solid circles,— triangles,— and 
squares— are the experimental data. 



NN hopping), U = 6.6, J z = -0.1, and Vi = 1.55. The 
adopted values for t', t" and U are about the same as 
those used in the literaturei&iil With the recently devel- 
oped algorithms for dealing with the summation over the 
Matsubara frequency and the Fourier transform of the ef- 
fective interaction with long-wavelength singularity^ we 
have numerically obtained the self-consistent solution to 
the Green's function. In this section, we present the nu- 
merical results for the physical quantities. 



A. Normal state 

In Fig. the numerical results (solid and dashed 
curves) are presented for Tn as functions of x. For 
comparison, the experimental data for NCC QVi 3 are de- 
picted by scaling the maximum T/v = 255 K to 0.0822£ 
(which means t w 270 meV). The present result for Tjv 
is about one order of magnitude smaller than that from 
the MFTiiii^ Even though the FLEX approach without 
consideration of the impurity scattering [vi = 0) gives 
the overall lower Tn, its non-monotonic behavior does 
not reflect the feature of the experimental data. In con- 
trast, the calculation with impurity strength Vi/t = 1.55 
fits very well the experimental data. In our calculation, 
the overall magnitude of Tn is controlled by U and J z ; 
larger U and J z lead to higher Tn- With increasing Vi 
from zero, the curve T/y(x) is suppressed gradually from 
the non-monotonic result toward a monotonic decreasing 
(as a function of x) one. 

In order to understand how the impurities suppress 
the AF transition, we firstly analyze the criterion 1 + 
Ux(Q) — in the absence of impurities. In this case 
r = 1, the spin susceptibility \ is simply determined 
by X(q,iu> n ). We note that at q = Q, a large con- 
tribution to x(Q) stems from the small regions around 
the hot spots kh which are the cross points of £g = 
and ££ + q = as shown in Fig. ^| For the free par- 





FIG. 10: Sketch of the hot spots in the first quadrant of the 
Brillouin zone. 



tides, one gets G^(k h ,iuj n )G^(kh + Q,iu> n ) = ~ w n 2 
since the energy £g = £,j: h+ Q = 0. In the present case 
the electrons interact via spin fluctuations, V s (q) has 
a divergent peak at q = {Q, 0) in the limit T — > Tjv 
from above. Considering only the predominant contri- 
bution from long- wavelength spin fluctuations, we then 
have approximately S(fc) ~ r ) 2 G{k + Q,iuj n ) with 7 2 = 
~3 2o V s {q)/'ZNf3 where the (^-summation runs over a 

small region around Q. By this approximation, G can 
be obtained analytically. Especially, at the hot spots, we 
have G(kh,ik> n ) ~ 2/ia) n (l + + 47 2 /o; 2 ), which im- 
plies a reduction in the magnitude of G. Even at large 
doping where hot spots vanish, the magnitude of G can 
be reduced by the spin fluctuations as well. Thus, the 
magnitude of X and thereby |x(Q)| are reduced. Since 
|x(Q)| increases with decreasing T, one therefore obtains 
a lower T/v than that from the MFT. In the presence of 
impurities, the magnitude of x(Q) is further decreased 
by r (|r| < 1), which leads to a further lowering of T/v. 

Here we also wish to give the reason why the FLEX 
result for T/v at Vi = is a non-monotonic function of x 
in the absence of impurity scattering. Again, we consider 
the criterion X+Ux{Q) — that determines T/v. As men- 
tioned above, a large contribution to x(Q) stems from 
the small regions around the hot spots kh- Let us look 
at the two hot spots in the first quadrant of the Brillouin 
zone. With increasing x from zero, these two hot spots 
move toward each other. Meanwhile, between the two 
hot spots, the two curves £g = and — approach 
each other too. Therefore, the inter-region enclosed by 
the curves between the two hot spots begins to give con- 
siderable contribution to x(Q)- This additional contribu- 
tion enhances the magnitude of x(Q)i ancl thereby leads 
to increasing of T/v with x. On the other hand, at large 
doping limit, there arc no hot spots. With increasing 
doping, the two curves depart from each other, resulting 
in a decreasing of T/v. Therefore, there is a maximum in 
Tjv at certain x within the FLEX calculation. With the 
impurity scattering, however, the density N% increases as 
the same as x, reducing the magnitude of the vertex func- 
tion r. Consequently, at small x, T/v cannot increase as 



FIG. 11: (Color online) Occupied spectral density around the 
chemical potential at temperatures T/Tv = 1, 1.18, and 1.84, 
at doping concentration x = 0.05. 

much as that predicted by FLEX. In the present case, 
with Vi — 1.55, Tjv monotonically decreases with x. 

It is expected that spin fluctuation effect may be ob- 
servable above but close to T/v- For example, the effect 
can be clearly reflected in the spectral density which is 
defined as 

A(k,E) = --lmG(k,E + i0+). (40) 

The variation of the spectral density with temperature 
may be detected by the ARPES experiment. By use 
of the Pade approximation^ the analytical continuation 
for the Green's function and thereby the function A(k, E) 
can be obtained. Fig. ^2 presents the occupied spectral 
density A<(k, E)) = f(E)A(k, E) (with / the Fermi dis- 
tribution function) around the chemical potential inte- 
grated within the energy window of — 0.02£ < E < 0.02t 
at T/T N = 1, 1.18 and 1.84, and at x = 0.05. At this 
doping, far below T/v, the spectral intensity is expected 
to mainly concentrate around the points (jr, 0) and (0, w). 
However, although there is no AF order, the Fermi sur- 
face splits into three parts at T above but close to T/v. 
The intensity loses at the hot spots, which is due to the 
formation of the pseudogap caused by spin fluctuations. 
When T above T/v continuously increases, the spin fluc- 
tuations become weak and the pseudogap diminishes at 
the hot spots. As a result, the Fermi surface becomes a 
continuous curve as seen in the case T/T/v = 1.84. 

The pseudogap can be also reflected in the den- 
sity of states (DOS) which is given by p(E) = 
(l/N)J2k A (k,E). Shown in Fig. UUare the results for 
DOS at x = 0.05 and various temperatures. As the tem- 
perature decreases from high value to T/v, DOS changes 
gradually. Especially, around E — 0, a visible depres- 
sion of DOS due to the pseudogap formation is seen as 
T — > T/v. The onset temperature T* for the depression 
at this doping is about 1.7T/v. 

In Fig. 1131 T* is shown as a function of doping con- 
centration, and is compared with the experimental data 
extracted from measurements on optical conductivity^ 
Clearly, our calculated T* is consistent with the experi- 
mental data. Similar result on T* has also been obtained 
by Kyung et alM> using a different approach. 

So far we showed the effects of spin fluctuations and 
impurity scattering on the properties of the normal state. 
In the next subsection, we present the results in the AF 
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FIG. 12: (Color online) Density ol states at x = 0.05 and 
T/Tjv = 1, 1.18, 1.48, 1.84. 
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FIG. 13: Pseudogap onset temperatures T* as a function of 
doping concentration x. The present result (theory) is com- 
pared with the experimental data of Ref. I27I 



state. 



B. AF phase 

We here begin with the spectral density in the AF state 
and compare it with the experimental data. Exhibited 




FIG. 14: (Color online) Occupied spectral density around the 
chemical potential at T/t = 0.003 and x — 0.1. 
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FIG. 15: (Color online) Occupied spectral density as function 
of momentum and energy at x — 0.1 and T/t — 0.003. The 
abscissa in (a)-(f) represent the momentum axis shown in Fig. 

ED 



in Fig. El is the occupied spectral density around the 
chemical potential at x = 0.1 and T/t = 0.003. It shows 
that the Fermi surface is near the two points (7r, 0) and 
(0, 7r) where the spectral intensity is strong. This kind of 
Fermi surface has been observed in ARPES experiment 
by Armitage et al~ Later, Matsui et al~ obtained the 
momentum-energy-resolved result for the occupied spec- 
tral density by ARPES experiments. In order to compare 
our result with the corresponding experimental data, we 
show the occupied spectral density in Fig. El as functions 
of momentum and energy. The momentum axes in Figs. 
I15f a)-(f) are along the lines a-f in Fig. El Figure H^T a) 
shows that the chemical potential cuts the upper band. 
The intensity in the region in the upper band close to the 
Fermi surface is very strong, illustrating the large weight 
of the quasiparticle there. On the other hand, the lower 
band is largely broadened, as shown by the wide area of 
weak intensity. From Figs. I15f a) to ll5f f). the intensity 
distribution shifts upward gradually, and the lower band 
starts to show stronger intensity. In Fig. I15f c). the inten- 
sity is nearly vanishing at the chemical potential. This 
corresponds to the case of line c in Fig. 1 141 where the line 
nearly misses to cross the Fermi surface. In Fig. I15f d) 
there is no intensity at the chemical potential since lind 
d as shown in Fig. ElP asses approximately the hot spot. 
Comparing these figures, we find that the gap between 
the chemical potential and lower band increases gradually 
from Fig. I15f a) to Fig. I15f where the largest gap shows 
up. In Fig. ITBT f). it is seen that the lower band is close 
to, but still away from the chemical potential. It should 
be pointed out that, at larger doping, the lower band 
will shift up and cross the chemical potential around the 
point (7r/2,7r/2). Correspondingly, the spectral intensity 
near (tt/2,tt/2) in Fig. 1141 will become strong enough 
and finally the several pieces of the Fermi surface will 
connect together. The present results are in fairly good 
agreement with the experimental observations!^ 

In Fig. El we exhibit the spectral density along the 
high-symmetry lines in the Brillouin zone within a large 
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FIG. 16: (Color online) Spectral intensity as function of en- 
ergy along the high-symmetry lines in the Brillouin zone at 
x = 0.1 and T/t = 0.003. 
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FIG. 17: (Color online) Density of states of AF phase at 
various doping concentrations at T/t = 0.003. 



energy scale to show the quasiparticle band structure at 
T/t = 0.003 and x = 0.1. In contrast to the mean field 
result, Fig. 1 161 shows no clear dispersion relation between 
the energy and momentum, but considerable broadening 
of spectral density. Only close to the chemical poten- 
tial the intensity shows sharp peak, implying that the 
quasiparticles have long life-time at the chemical poten- 
tial. At (7r/2,7r/2), the gap clearly shows up. The upper 
band is very flat at (jr, 0), where it reaches its minimum. 
At low doping, the chemical potential is close to the bot- 
tom of the upper band. Even at zero doping, due to the 
many-body effect, the upper band is already occupied 
with quasiparticles; the chemical potential is at nearly 
the bottom of the upper band. From such a band struc- 
ture, we can expect an unusual DOS. Shown in Fig. 1171 
is the DOS of low temperature AF states at several dop- 
ing concentrations. Clearly, a central peak appears in 
the DOS at zero energy at low doping, which originates 
from the fact that the chemical potential cuts the flat 
region near the point (tt, 0) where the width of energy 
level is very narrow. With increasing the doping concen- 
tration the Fermi surface departures from the flat region, 
the peak in DOS at large doping is less pronounced. At 
x = 0.125, the peak diminishes as shown in Fig. 1171 
and it is expected to disappear at higher doping concen- 
tration. At zero doping, since the Fermi surface is very 
small, the zero-energy peak almost diminishes too. It 
should be noted that since the direct gap as appeared in 
the spectral density varies with the momentum, we ob- 
tain a structured gap in the DOS that is a superposition 
of the spectral density. We here call the distance between 
the humps below and above the chemical potential in the 
DOS as the gap. In Fig. 17, the gap appears between 
E ~ — 0.5t and ~ 0.3t. Especially, the bottom of the 
upper Hubbard band falls into the DOS gap consistent 
with the situation shown in Fig. 1161 This is different 
from the MFT by which one obtains a DOS with a clear 
gap between the upper and lower bands, but not a peak 
at the chemical potential. 

What is the physics of the structured gap with a cen- 



tral peak at zero energy in the DOS? In the system, there 
can be two kinds of elementary excitations: the single 
quasiparticles and the spin waves due to the transverse- 
spin fluctuations. The energy for the latter is within a 
range < fl < fl c . Now, consider a quasiparticle with 
energy E > above the chemical potential. The quasi- 
particle can excite a spin wave by transferring energy £1 
to the latter and transiting itself to a low energy E — 
state. Such an excitation can happen in the energy range 
< fi < min(O c , E). Thus, the lifetime of the quasipar- 
ticle is finite, and the DOS at E is suppressed. On the 
other hand, for the particle at very low energy, the en- 
ergy range of the excited spin waves is very restricted; 
the integrated effect of the spin-wave excitations on the 
particle is very small. It means that the quasiparticles 
at very low energy have long lifetime. Therefore, the 
spin-wave excitations result in a dip in the DOS of the 
single particles above the chemical potential, but leave 
a peak at zero energy. The similar mechanism of the 
peak-dip-hump structure in the DOS of the HTS due to 
a dispersionless collective mode has been explained by 
Shen and Schriefferi2& It should be pointed out that be- 
sides the well-defined spin waves at low-energy regime, 
the spin fluctuations at higher energy can modify the 
DOS as well; the effects from both of them are included 
in the present calculation. The dip can be understood 
as pseudogap. The gap between the lower band and the 
chemical potential which is at about the bottom of the 
upper band corresponds to the mean field AF gap caused 
by the AF order. In the present case, it is rounded by 
the spin fluctuations. 

The present zero-energy peak in DOS may correspond 
to the in-gap states investigated by Kusunose and Rice. 9 
The origin of a in-gap state was explained as the admix- 
ture of the quasihole mean filed state in the lower band 
with the hole state accompanied by particle-hole exci- 
tations due to spin fluctuations in the upper band. By 
adopting the self-consistent Born approximation with the 
RPA spin-spin interaction and the MFT order parame- 
ter, they obtained the in-gap states at energies about 
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FIG. 18: (Color online) Distribution function n(k) in the first 
quadrant of the Brillouin zone at x = 0.1 and T/t = 0.003. 

—t away from the chemical potential. In the present 
case, however, the order parameter, the spin fluctuations 
and the chemical potential adjust themselves all in a self- 
consistent way. As a result, the peak appears always at 
the chemical potential. 

In Fig. 1181 we depict the momentum-space distribution 
function n(k) (at x — 0.1) defined by 

n(k) = — ^2 G{k, iuj n ) exp(iw n r)) (41) 

n 

where 7/ — > + is an infinitesimal small positive constant. 
Even at very low temperature T/t — 0.003, the distribu- 
tion has a small but finite width without abrupt discon- 
tinuity at the Fermi surface, and thus is non-Fermi-liquid 
like. In addition, within the Fermi area, the states are 
not fully occupied in contrast to the case in a free electron 
system. Correspondingly, outside the Fermi area n(k) is 
finite; the minimum of n(k) is larger than 0.1. The dis- 
tribution n(k) shown in Fig. 1161 is a typical result within 
the doping range < x < 0.17. 

Figure ITTjl exhibits the staggered magnetization m as a 
function of x at very low temperature T w §K, and the 
comparison with the experimental dataA^S The exper- 
imental value for m(0) at Cu 2+ is about 0.5 in NCCO 
and ~ 0.4 in PCCO, 4 which are close to our theoreti- 
cal result 0.48. The overall behavior of m versus x ob- 
tained by the present calculation is in very good agree- 
ment with experiment. Though the dynamical- MFT cal- 
culation by Sencchal et alm& yields a reasonable behavior 
for m(x)/m(0), their result m(0) w 0.7 is too large. 



IV. CONCLUSION AND REMARK 

In conclusion, we have presented the formulation of 
antifcrromagnetism in electron-doped cuprate supercon- 
ductors on the basis of the Hubbard model using the 
fluctuation-exchange approach. A self-consistent treat- 
ment of the spin fluctuations in combination with the 
impurity scattering due to the disordered dopant-atoms 
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FIG. 19: (Color online) The calculated staggered magnetiza- 
tion m (solid line) as a function of doping concentration x at 
T fa OK. The experimental data are also shown for compari- 
son [solid and hollow circles^- solid triangles—]. 



can reasonably account for a number of experimental re- 
sults including the AF transition temperature Tjy, the 
onset temperature of pseudogap T* , the spectral density, 
and the staggered magnetization m. The single parti- 
cle distribution function has also been calculated, and it 
differs significantly from that of the Fermi liquid. Also, 
we have found the existence of zero-energy peak in the 
density of states of the AF phase. The peak is sharp in 
the low doping regime, and disappears near the optimal 
doping. This result is not inconsistent with the ARPES 
experiments, but still needs to be carefully checked by 
the tunneling experiments. 

Although the present theory provides reasonable ex- 
planations for a number of experiments, our calculated 
spin- wave velocity c, however, is smaller than the ex- 
perimental data. With the hopping parameter t ss 270 
meV, the experimental results for c are 0.74ia and 0.94ta 
respectively for PCCO and NCC02S (a: the lattice con- 
stant of CuO plane), while our theoretical prediction is 
only about 0.33ta. The quantity c is related to second or- 
der derivative of the transverse spin susceptibilities with 
respect to the momentum (Px^iQ + q)/dq 2 \ q=0 . In the 
present approach, we are probably able to reasonably 
treat the overall behavior of the two-particle propagator 
Xui,(<?)> but not its second order derivative with respect 
to q. Further study is needed in this point. 
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APPENDIX 

In this Appendix, we deal with the average of V/; 1 (q) 
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over —it < q z < ir. As illustrated in Section II, we ne- 
glect the weak g z -dependence in the numerator. Here, we 
consider only the case of Q z = 0. The result for Q z = it 
can be derived analogously. 

For Q m — 0, Xi2(l) = 0) the expression for Vi 1 (q) 
reduces to the similar form as Eq. (flfif) . Therefore, the 
average Vyi(q) is given by 



u 2 x\M 



For 17 m 7^ 0, the denominator of Vl±{q) that is the 
determinant D — det\\ + 'x t qV{q)\ can be written as D = 
a cos 2 q z — b cos q z + c with 

c = l + Ulx^+xU^+aiU/Jz) 2 . 

Notice that X12 (?) is a pure imaginary quantity. D never 
vanishes at Q m ^ 0. By denoting the average by 



the function g(q) is given by 
(A) d = b 2 - Aac> 0, 



5(9) = -7=[sgn(a;_) — - sgn(x+) 



where x± = (b ± \fd)/2a. 
(B) d = 0, 



5(9) 



a(.T 2 - l) 3 /2 



where x = b/2a. 
(C) d < 0, 



? (g)=-^=Im(z 2 -ir 1 / 2 



yA(9) = ^ 2 (i + a^/^ 2 )ff(g), 



where z = (\b\ + iy/\d\)/2 
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